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A detailed comparison is made between the topological structure of quenched QCD as revealed by 
the recently proposed over-improved stout-link smearing in conjunction with an improved gluonic 
definition of the topological density on one hand and a similar analysis made possible by the overlap- 
fermionic topological charge density both with and without variable ultraviolet cutoff Xcut- The 
matching is twofold, provided by fitting the density-density two-point functions on one hand and by 
a point-by-point fitting of the topological densities according to the two methods. We point out the 
similar cluster structure of the topological density for moderate smearing and 200 MeV < Xcut < 
^ ■ 600 MeV, respectively. We demonstrate the relation of the gluonic topological density for extensive 

lO ! smearing to the location of the overlap zero modes and the lowest overlap non-zero mode as found 

■ for the unsmeared configurations. 

PACS numbers: ll.15.Ha, 12.38.Gc, 12.38.Aw 

^ ■ 

o . 

oo I- INTRODUCTION 

o" 

Massless overlap fermions [1, 2] have provided us with a fermionic definition of topological charge [J, 

■ This offers the advantage that, by truncating to the overlap modes with |A| < Xcut^ the effect of 
rS I ultraviolet filtering can be studied [5] without changing the gauge field itself. A broad study of 

^ ■ different aspects of vacuum structure, so far for quenched QCD, has been published recently [6]. On 

the other hand, during the 80 's procedures of cooling or smearing of gauge fields were proposed to 
exhibit the infrared structure of gauge fields 0, • This has been considered with reservations because 
it is difficult to assess in which aspect the gauge field configuration could have changed under this 
procedure. The practitioners of cooling/smearing, in particular those who were focussing on vacuum 
structure in the form of extended, smooth (semiclassical) structures, have continued to improve their 
techniques: cooling with improved actions [9], restricted improved cooling [l^ etc). Over-improved 
cooling [HI has been applied in order to prevent instantons or other topological excitations from 
collapsing. 

More recently, the concept of cooling/smearing has been modified in another direction, to make 
it differentiable with respect to the unsmeared field. The so-called stout-link smearing [13] makes it 
possible to evaluate the "force" for the molecular dynamics in non-quenched update algorithms with 
a fermionic action involving some form of "fat" (smeared) links substituting the original gauge links. 
The concept of over- improvement has also recently been applied to stout-link smearing [14]. 
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In this paper, following Ref. [15] we want to characterize the subsequent mapping "link stout link" 
as a particular relaxation scheme, eventually leading to a finally classical configuration. One obvious 
way to discuss this process is to record the local distribution of topological charge density expressed 
by an improved [16] gluonic definition of field strength and topological density. In this paper we shall 
compare the emergent structures with the topological density provided by the overlap definition [6| 
with different levels of ultraviolet filtering. Surprisingly, the overlap topological charge density without 
filtering, that recently has been found to form lower-dimensional structures 0, [T^^ES^ corresponds to 
few steps of smearing. This comparison will lead us to a one-to-one mapping of the ultraviolet cutoff 
Xcut (mode truncation in the overlap picture) to the number of stout-link smearing iterations over 
a wide range of smearing iterations. In a similar spirit, the correspondence of APE smearing ^|, 
Laplacian filtering [20] and the topological density filtered according to another Ginsparg- Wilson Dirac 
operator [21, 22] has iDeen studied recently 0, [13]. 

Only a few iterations of stout-link smearing are necessary before structures become recognizable 
with the gluonic definition of the topological density, and these structures are surprisingly far from 
4D extended, sign-coherent lumps, such that the topological density compares well with the unfiltered 
overlap definition. We have two criteria to establish this relation between smearing and filtering. First, 
it is the behavior of the two-point correlation function of the topological density that emerges from 
the respective definition. The second is the actual site-by-site difference of the topological density 
profile over a set of lattice configurations. The quality of the latter coincidence is surprisingly good, 
which supports the reliability of both methods to explore the vacuum structure. 

Finally, however, this relation becomes loose because stout-link smearing turns the configurations 
into piecewise classical fields which apparently resemble instantons and anti-instantons. Although 
the coherence among the lowest overlap modes guarantees a relatively simple picture of the fermionic 
topological density and the ultraviolet filtered gluonic field strength with low cutoff Xcut [6], there is 
no argument as to why the overlap picture should be in correspondence to a (linkwise) classical lattice 
configuration. 

It is intriguing to see that the instanton-like structure that is revealed in this late stage of smearing 
corresponds to the overlap zero mode(s) and the lowest pair(s) of overlap non-zeromodes obtained for 
the respective unsmeared (equilibrium) configuration. A similar observation has already been made 
by Negele et al. [25]. In Ref. [18] some of us have presented a cluster analysis of individual eigenmodes 
of the overlap Dirac operator. It turns out that the moderate number of clusters that the zero and 
first non-zero modes consist of (at a level of scalar density below the peak values) are pointing towards 
the positions where instantons and anti-instantons appear later. 

Our paper is organized as follows: in Section[IIlwe explain briefiy the overlap operator and the stout- 
link smearing procedure, in Section IIIII we describe the matching between smearing and ultraviolet 
filtering, according to the two-point correlator and according to a global fitting of the profile of charge. 
In Section IIVI we try to relate the clusters that both definitions exhibit to each other. In the final 
stadium of smearing we shall confront the emergent semiclassical lumps with the lowest eigenmodes 
(zero mode and lowest non-zero mode) of the original configurations. In Section [Vl we draw conclusions 
and give an outlook. 

II. TECHNICAL DETAILS 
A. Configurations and the overlap definition of topological density 

The configurations underlying this comparison stem from an extended investigation published in 
Ref. 0]. The configurations are taken from an ensemble of 16^ x 32 quenched lattices generated with 
the tadpole improved Liischer-Weisz action at = 8.45. 

The (massless) overlap Dirac operator is constructed for the Wilson-Dirac input kernel Dw = 
M — p/a, M being the massless Wilson-Dirac operator with r = 1, and p = 1.4. The corresponding 
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solution of the Ginsparg- Wilson relation reads as follows 

^(0) = ^ (l + Dw/^D\^Dw^ = a ^ i^sgii{Hw)) , (1) 

with Hw = jbDw Circa 150 overlap eigenmodes have been obtained per configuration. They have 
been used to construct ultraviolet smeared topological densities according to a cut-off Xcut = 200 MeV, 
400 MeV and 635 MeV. For half of the subset of 10 configurations particularly considered in the 
present study we have also calculated for Refs. [l^ the overlap topological density without mode 
truncation, the "all-scale" topological density. 

The spectrum eventually consists of some zero modes, in addition to pairs of non-zero modes of 
globally vanishing chirality. The topological density can be formally obtained from the trace of the 
overlap Dirac operator [4] 

q{x) = -tr [75 (1 - ^1^(0; x, x))] . (2) 

Using the spectral representation of the overlap Dirac operator, a family of ultraviolet filtered topo- 
logical charge densities labelled by qx^^^ {x) can be obtained [H, [26| , 

QXout{x) = - Yl (^-^)^x(^h^^x{x). (3) 

The topological charge of each configuration fulfills the Atiyah- Singer index theorem [3] 

Q = n_ - n+ . (4) 

Note that the zero modes of any given configuration carry the same chirality. 

In this paper we shall point out that the above-mentioned family of densities is well represented by 
the gluonic topological density after an appropriate number of iterations of stout-link smearing. 



B. Over- improved stout-link smearing 



In Refs. P^, ll5| two of us have motivated the stout-link smearing with respect to a specific over- 
improved action. The concept of over-improved action arose from the observation that the Wilson 
one-plaquette action does not guarantee the stability of instantons [2?]. Even with the Symanzik 
improved action including 1x1 plaquettes and 1x2 rectangular loops, for a classical instanton 
solution of size pinst placed on the lattice, the action is obtained as 



-yinst 



n2 



17 



210 \ pinst 



(5) 



Following Garcia Perez et al. [11], Ref. [14] proposes to modify the relative contribution of plaquette 
and rectangular terms, parametrized by a parameter e. For the Symanzik action this reads as follows 



2e 



(1 - P^..{x)) 



1 



(1 - R^^ix) + 1 - R^,^{x)) 



(6) 



interpolating between Wilson (e 
plaquettes are denoted as 



3 ^ 12 
= 1) and Symanzik action (e = 0). As usual, = 6/^^. The 
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and the traces of the rectangular loops as 
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With this action the instanton action behaves as 



o 2 



^ e f a \_ Ue-17 



5 \pinstj 210 \pinst 



(7) 



In Ref. [15] a value e = —0.25 was found to provide an instanton action almost independent of the 
radius for pinst > 1.5a. 

The stout-link iteration solves the equation of motion (without momentum) 



U"J,' - C/rr = exp ( -p,^ ^ ) C/°|f (8) 



in parallel, with a "time step" psm- No projection onto the group 5/7(3) is required in distinction to 
APE smearing. The derivative is the traceless antihermitean part of a sum of untraced loops 

S S (c^ 5 26 1 6 

-^fj^ = Yl (^^^ ^^^P^ touching U) + -j^ Yl (1^2 + 2x1 loops touching U)^^^,.^^^^ ^^^^^^^^^ 

with Ux^ in the leftmost position. This method has been discussed in Ref. [lH with respect to the two- 
point correlator of the topological density. The differences between quenched and dynamical ensembles 
of configurations with respect to the distribution of topological density have been demonstrated. 
The instanton-like features are seen to become stronger with more and more iterations of (|8|). The 
resemblance to instantons was stated [15] by a close correlation - maximum by maximum of the 
modulus of the topological density - between the "instanton size", determined by the curvature of 
\q{x)\ in the maximum position xq of the cluster, and the value of the density |<7(xo)|. 

The local lattice operator [16] to represent the topological density for the stout-link smeared con- 
figurations is based on a highly improved field-strength tensor, 

F;r(a:) = hClYHx) + k^Cf^Kx) + k,C''^^\x) + k,C^^f\x) + k,Cff\x) (9) 

with 



fcl = 


19 

~9 ~ 


55/C5 , 


k2 = 


1 

36 ~ 


16^5 , 


ks = 


32 

~45 
1 

15 ~ 


+ 64/C. 


= 


6/C5 , 



(10) 



and with 



clr\^) = I {wi,:'"^\x) + w;,T"H^)) (n) 

being a symmetrized "clover sum" of {nx p^mx 0) Wilson loops around the site x. A 3-loop improved 
field strength tensor can be achieved chosing = 1/90, such that /cs = /c4 = 0. The topological charge 
density is then represented in the form 

2 

qsmix) = ^Tr (f^uF^^) ■ (12) 
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n,^ Q = -8 Q = Q Q = Q = -7 Q = Q 



1 


0.9954 


-0.0001 


1.0054 


0.6442 


0.0000 


2 


0.9634 


-0.0001 


-21.5369 


0.7987 


-0.0001 


3 


0.9494 


-0.0001 


3.6497 


0.8990 


-0.0001 


4 


0.9585 


-0.0002 


1.8690 


0.9485 


-0.0001 


5 


0.9760 


-0.0004 


1.4537 


0.9728 


-0.0001 


6 


0.9903 


-0.0008 


1.2736 


0.9843 


-0.0002 


7 


0.9982 


-0.0015 


1.1725 


0.9894 


-0.0004 


8 


1.0010 


-0.0030 


1.1105 


0.9915 


-0.0006 



TABLE I: Table of Zsm values for the initial five configurations. The normalization procedure is only valid 
for Q 7^ and works best for large |Q|. The fluctuating values for the configuration with Q = — 1 occuring at 
T^sw 2 are due to Qsm being approximately at this point. 



III. MATCHING STOUT-LINK SMEARING TO OVERLAP FILTERING 
A. Matching the two-point correlator 

The two-point correlators for stout-link smearing and overlap filtering are matched using a mini- 
mization of the sum of the absolute difference between the correlators. That is, we compute 

min I (9(x)g(0)),„ - (g(x)9(0))A.„, 1^ (13) 

as a function of the number of smearing sweeps, Usw^ for fixed Xcut- 

The topological charge of each overlap filtered configuration is integer valued because they satisfy 
the Atiyah-Singer index theorem (|4]). For stout-link smearing it can take up to 10 sweeps of smearing 
to achieve an integer charge on these 16^ x 32 lattices. In order to compare overlap filtering and 
stout-link smearing for Usw < 10 in a fair way a non-perturbative normalization is applied. Given 
that we know the topological charge Q from the overlap configurations, we calculate, for each number 
of smearing sweeps Usw^ a normalization factor Zsm via 

Q = Zsm^Qsmix) , (14) 

X 

and Qsmix) is then normalized through qsm{x) Zg^n qs^{x). This ensures that Qsm = Q- We also 
experimented with an alternate normalization where we matched the absolute values of the topological 
charge density, however this proved less fruitful. 

Typical values of Zsm are provided in Table [D For obvious reasons, it is only possible to extract 
a Zsm factor for Q 7^ 0, and the procedure works best for Q far from zero. Consequently, the best 
results are found for the configurations where Q = —7 and —8, where the Zsm values rapidly approach 
1. This is also the case for the configuration with Q = — 1, however Zsm fluctuates at Usw = 2, which 
is due to Qsm being approximately at this point. 

In Table [TT] we present the best matches for the two-point correlators between the filtered overlap 
densities and the stout-link smeared gluonic densities. Configurations have been used for which the 
unfiltered overlap topological density (full density) has been measured. The correlators of the fermionic 
topological density including all known modes, with Xcut = 634 MeV, with Xcut = 400 MeV and with 
Xcut = 200 MeV are matched against the smeared gluonic correlators. 

We show in Fig. 1 the two-point correlator of the topological density for a single configuration 
represented by the unfiltered fermionic topological density (full density) of Eq. ([2|) of the equilibrium 
configuration on one hand and for the gluonic definition of Eq. (fT2]) after 5 smearing steps. It is 
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= -8 


= 


Q = -1 


Q = -7 


Q = 


full density 


005 


005 


005 — 


005 005 


005 004 


005 — 


all known modes 


039 


039 


043 — 


036 036 


038 037 


034 — 


634 MeV 


052 


052 


064 — 


048 048 


055 055 


048 — 


400 MeV 


127 


127 


183 — 


109 108 


139 138 


115 — 


200 MeV 


248 


247 


300 — 


232 232 


300 300 


271 — 



TABLE II: The best matches for the two-point correlators of five different configurations as determined by 
Eq. (p^ . For each configuration, the two columns give the number of smearing sweeps that correspond to 
the best match for a given level of UV filtering. In each case the left column gives the best match for the 
unnormalised gluonic density and the right column gives the best match for the gluonic density normalised 
with Zsm (see text). The right columns for the Q = configurations are absent because the normalization 
procedure fails in these cases. 
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FIG. 1: The two-point function of the fermionic topological density for one configuration shown without mode 
truncation compared with the bosonic definition after 5 steps of smearing. 

remarkable that the two correlators follow each other's fluctuations at larger distance. There is some 
difference in normalization of the negative peak. This configuration has Qferm = —8, and the gluonic 
definition gives Qsm = —8.19. In this case the "non-perturbative" renormalization would even slightly 
increase the difference between the curves as the three curves in Fig. 2 show. 

In Fig. 3 the density-density two-point correlator for the same configuration is compared for an 
ultraviolet cutoff Xcut — 634 MeV for the fermionic, overlap definition and the gluonic definition after 
52 smearing steps. The filtered correlators match perfectly. The third curve shows the correlator with 
the gluonic definition after 45 smearing steps. This refers to the case of an optimal point-by-point 
matching of the densities, as discussed in the next section. The corresponding correlator is somewhat 
higher and steeper because it uses slightly less smearing iterations. 

Fig. 4 shows the same configuration for a lower ultraviolet cutoff \cut = 400 MeV. After 127 
smearing steps the correlator for the gluonic definition fits the correlation function perfectly. Somewhat 
less smearing steps (103), again optimally fitting the fermionic topological density point by point, tend 
to overestimate the correlator. 

All smeared correlators discussed above were generated using a smearing parameter of psm = 0.06. 
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FIG. 2: The same as in Fig. 1 with the effect of the renormahzation factor Zsm = Q ferm/Qsm also shown. 

(Qferm — 8, Q sm — 8.19) 
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FIG. 3: The two-point function of the fermionic topological density with an UV cutoff at Xcut — 634 MeV 
compared with the bosonic definition after 52 steps of smearing when the correlator is fitted best. Smearing 
after 45 steps leads to the best global matching of the density qx^^^ with Xcut — 634 MeV, but the correlator 
is steeper. 



By using a different value for psm it is possible that a different number of smearing sweeps will 
provide the best fit. This is because the amount of smearing applied to a gauge field is proportional 



to ps 



By varying psm as well as the number of sweeps one has greater fine-grained control 



over the matching. Using a variable Psm, but holding Usw = 5 fixed, and applying this to the same 
configuration considered previously we find that psm = 0.055 provides the best match for the unfiltered 
topological density. 
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FIG. 4: The same as Fig. 3 with an UV cutoff at Xcut = 400 MeV fitted best by the bosonic definition after 
127 steps of smearing. Smearing after 103 steps leads to the best global matching of the densities, but the 
correlator is steeper. 
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FIG. 5: The same as Figs. 3 and 4, now with an UV cutoff at Xcut = 200 MeV fitted best by 248 sweeps of 
smearing. 260 sweeps of smearing gives the best match for the topological charge densities. 



A comparison of some different values for for the smearing parameter is provided in Fig. [6l As 
expected, increasing psm results in a suppression of the negativity of the two-point correlator. Al- 
though psm = 0.055 provides the best match through a minimization of the sum of the differences in 
the two-point functions, psm = 0.06 gives the best match for the x-intercept. The two-point function 
for this larger psm also appears to have a similar shape to the two-point function of the full fermionic 
topological density. 
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FIG. 6: A comparison of the smeared two-point function for three different values of the smearing parameter, 
Psm- The configuration shown is the same as that used in the previous figures. A best match as dictated by 
Eq. (|13|) is given by psm — 0.055, however psm — 0.06 gives the best match for the x-intercept. 





Q = -8 


Q = 


Q = -1 


0= -7 


Q = 


full density 


300 005 


300 — 


300 006 


300 005 


300 — 


all known modes 


037 036 


037 — 


035 035 


035 035 


033 — 


634 MeV 


045 045 


048 — 


044 044 


048 048 


044 — 


400 MeV 


103 103 


098 — 


103 103 


087 087 


093 — 


200 MeV 


261 260 


300 — 


236 236 


192 191 


187 — 



TABLE III: Best matches for q[x) between the filtered overlap densities and smeared gluonic densities. The 
left columns contain the best matches with the unnormalised gluonic densities, and the right columns contain 
the best matches when including the Zsm factor. We see the importance of the Zsm normalization factor when 
attempting the match the full fermionic density. The right columns are absent for the Q = configurations 
because the normalization method does not work for this Q. 



B. Matching the fermionic and gluonic topological density point by point 

We now aim to match the filtered fermionic densities with the smeared gluonic densities through 
a point-by-point matching of the respective topological charge densities qx^^t {x) and qsm {x) . Given 
some filtered fermionic topological charge density we compare it against the gluonic density Qsmix) 
after some number of sweeps by calculating the absolute value of the difference between the respective 
q at each x. To find the best match we compute the minimum of the sum of the differences, 

(^\Qsm{x) - qx,^,{x)\^ , (15) 

as a function of the number of smearing sweeps, risw^ for fixed Xcut- The non-perturbative normaliza- 
tion of Eq. ([llj) will also be applied. The best matches are presented in Table Hill For a small number 
of sweeps the importance of the Zgrn normalization is apparent. Visualizations of some configurations 
are shown later in Section HVl 
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C. Towards the no-smearing limit 



Increased smearing leads to flattening of the Euclidean two-point function. It follows that less sweeps 
of smearing leads to an increasing negative dip in the correlator, and we now study the behaviour 
of the two-point function in the limit risw 0. After only one or two sweeps of smearing there is a 
non-trivial renormalization that must be applied to the topological charge density. From Table [J we 
see that after three smearing iterations the Zsrn factors for the Q = —7 and Q = —8 configurations 
are in relatively good agreement. We can therefore study the correlator for Ugw ^ 3, averaging over 
these two configurations. 

Such a comparison leads to a series of correlators which are displayed in Fig. H The x-intercept 
also moves further towards zero. 
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FIG. 7: The two-point function as a function of smearing sweeps. The correlator has been averaged over the 
two configurations with Q = — 7 and Q = — 8. Both charge densities were normalised by an average of the 
Zsm values for the two configurations. As the number of sweeps is decreased the negative dip increases in 
magnitude and the x-intercept moves further toward zero. 



IV. COMPARISON OF TOPOLOGICAL CLUSTERS 



A. Clusters of both topological densities compared for weak stout-link smearing 

Using the matching of Eq. ([15]) we are able to directly compare the overlap topological charge density 
with some level of UV filtering to a given number of stout-link smearing sweeps. In an early stadium 
of stout-link smearing the topological density does not yet show classical, instanton-like features. 
What is meant by "instanton-like features" and how they gradually emerge from stout-link smearing 
is illustrated for a configuration with Q = — 1 in Fig. [8l The solid lines represent the relation 

\<lM\ = , (16) 

^ Pinst 

typical for the (anti)instanton solution, between the gluonic topological charge density in the maxima 
xo of the modulus of the density \qsm{x)\ and the "instanton radii" pinst obtained from a fit of the 
curvature of the action density in the points neighbouring xq. The upper plot has to be considered 
with a grain of salt because a closer look at the maxima reveals a substantial lack of isotropy of the 
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FIG. 8: Three scatterplots showing the increasing instanton-Uke correlation between the gluonic topological 
charge density in the maxima xo of its modulus \qsm{x)\ and the "instanton radii" pinst (see text). The solid 
lines represent the (anti) instanton- like relation between the two cluster parameters. The upper plot shows a 
huge number of maxima after 5 smearing steps, without any relation between density and size. The middle 
and bottom plots show a decreasing number of maxima and an inreasing accuracy of the instanton- like relation 
after 40 and 200 smearing steps, respectively. 
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FIG. 9: The gluonic topological charge density very close to the maxima in some timeslice of the same 
configuration as analyzed in Fig. [8] after 5 (left) and 40 stout-link smearing iterations (right). In color online: 
negative density blue/green, positive density red/yellow. In grey-scale: negative density dark, positive density 
fight. 




FIG. 10: The fermionic topological charge density of a Q = configuration with Xcut = 634 MeV (left) 
compared with 48 sweeps of over-improved stout-link smearing (right). In color online: negative density 
blue/green, positive density red/yellow. In grey-scale: negative density dark, positive density light. 



peaks of topological density which is however implicitely assumed in the fits of Pinst- The marked 
difference in this respect between 5 and 40 smearing steps is visible in Fig. [9] showing a certain timeslice 
of the same configuration as in Fig. El 

For a moderate amount of smearing and filtering, respectively, the two topological density definitions 
are faithfully exhibiting the outstanding clusterization of charge, provided the cut-off Xcut and the 
number of stout-link smearing steps are optimally tuned to each other. This is exemplified by Fig. [10] 
which shows the same time-slice of a Q = configuration, on the left-hand side portrayed by the 
overlap-fermionic topological density with Xcut = 634 MeV and on the right-hand side by the gluonic 
topological density after 48 stout-link smearing steps. 

In Fig. [11] we show the fermionic density of the Q = —8 configuration, with Xcut = 400 MeV and 
the gluonic density after 103 stout-link smearing steps. In this stadium of smearing, the bias in favor 
of classical lumps is already visible. 

Besides the similarity between the two methods, the tendency towards classical lumps becomes 
obvious in Fig. [12] This figure shows the Q = —7 configuration with a cutoff of Xcut = 200 MeV. 
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FIG. 11: The fermionic topological charge density of a Q = —8 configuration with Xcut = 400 MeV (left) 
compared with 103 sweeps of over- improved stout-link smearing (right). In color online: negative density 
blue/green, positive density red/yellow. In grey-scale: negative density dark, positive density light. 





FIG. 12: The fermionic topological charge density of a Q = —7 configuration with Xcut = 200 MeV (left) 
compared with 192 sweeps of over-improved stout-link smearing (right). In color online: negative density 
blue/green, positive density red/yellow. In grey-scale: negative density dark, positive density light. 



At this level of UV filtering the best match is provided by 193 sweeps of over- improved stout-link 
smearing. In this stadium of smearing, in particular because of the over-improved action built into 
the smearing, the minority positive charge has become stabilized compared with what the fermionic 
view shows. 



B. Zero modes and lowest non-zeromodes and the instanton content after long smearing 



For long smearing with essentially more than 100 smearing steps, the similarity of the topological 
densities slowly becomes less perfect. Still, the position of the gluonic topological lumps is not com- 
pletely arbitrary from the point of view of the original configuration. The zero modes (if Q 7^ 0) and 
the lowest non-zero modes contain a high predictive power over where these lumps will appear. 

For 200 smearing steps, this is illustrated in Fig. [13] by a Q = — 1 configuration. In the left panel 
the gluonic topological density is shown, in the right panel the scalar density of the zero mode. The 
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FIG. 13: The gluonic topological charge density after 200 smearing sweeps (left) and the scalar density of 
the zero mode (right) for a Q = — 1 configuration. One can see how the zero mode extends over only regions 
of negative charge. All extended sign-coherent objects seen are good local approximations to instantons or 
anti-instantons in the centre. In color online: negative density is blue/green, positive density is red/yellow 
(left), and scalar density is blue/green (right). In grey-scale: negative density dark, positive density light 
(left). 




FIG. 14: The gluonic topological charge density after 200 smearing sweeps (left) and the scalar density of the 
zero mode (right) for a Q = 1 configuration. In this case the high density regions of the zero mode are centred 
on lumps of positive charge. Again all objects are good approximations to classical instantons. In color online: 
negative density is blue/green, positive density is red/yellow (left), and scalar density is red/yellow (right). In 
grey-scale: negative density dark, positive density light (left). 



zero mode covers three distinct centers of topological charge of appropriate sign only in the selected 
time-slice. Two more examples of total charge Q = ±1 are shown in Figs. [T4l and [T5l One sees that 
the zero mode does not always cover all regions of appropriate charge. In other words, the gluonic 
version of the topological charge density for some clusters - even in a late stadium of smearing - may 
be built by non-zero modes. This leads us to revise the naive expectations according to which each 
zero mode would be residing on one lump of excess topological charge, for which a cluster charge of 
Qci = ±1 would be suggested by the instanton model. In fact, in Ref. [18] it has been demonstrated 
that the zero modes typically are simultaneously carried by a moderate number of centers below the 
peak density, before they start percolating throughout the whole lattice at an even lower level of the 
scalar density. Here we additionally learn that all theses centers are marked by lumps of topological 
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FIG. 15: The gluonic topological charge density after 200 smearing sweeps (left) and the scalar density of the 
zero mode (right) for another Q = — 1 configuration. The high density regions of the zero mode are again 
centred on lumps of negative charge. All topological objects are good approximations to classical instantons. 
In color onUne: negative density is blue/green, positive density is red/yellow (left), and scalar density is 
blue/green (right). In grey-scale: negative density dark, positive density light (left). 



charge of appropriate sign, however not necessarily all lumps are covered. 

For these same three configurations we now consider the distribution of the lowest non-zero modes. 
Starting with the first Q = —1 configuration, we present the scalar and pseudoscalar densities of the 
first non-zero mode in Fig. [161 One sees that the pseudoscalar density, according to its local chirality, 
highlights certain parts of the topological lumps with appropriate sign of charge, and leaves others 
(for other low- lying modes). 

The densities for the next two Q = 1 and Q = —1 configurations are shown in Figs. [T71 and [T8l 
In the instanton model, a whole band of almost-zero modes is generated by diagonalizing the Dirac 
operator in the field of a superposition of n+ instantons and n_ anti-instantons in the basis of linear 
combinations of the n+ + n_ zero modes corresponding to the case of infinite diluteness. Apart from 
the remaining \Q\ zero modes, the almost-zero modes are expected to bridge at least one pair of 
instanton and anti-instanton with the scalar density peaking on top of the topological charge lumps 
and the pseudoscalar density peaking with the appropriate sign there. Qualitatively, this is visible in 
Figs. [HIlTl and [181 

V. CONCLUSIONS AND OUTLOOK 

In this paper we have confronted the overlap-fermionic topological charge density and the improved 
gluonic topological charge at different levels of ultraviolet smoothing, realized in one case by a trunca- 
tion of the mode expansion at Xcut and in the other case by a certain number of smearing steps applied 
to the gauge field in order to wipe out ultraviolet fluctuations. These two views of getting the infrared 
topological structure of the gauge field correspond to each other. A similar result, however for APE 
smearing and a different improved Dirac operator, has been found in Ref. |23|,|2J]. In the present paper 
this has been confirmed for two other specific realizations of both methods, using the massless overlap 
Dirac operator Q on one hand and stout-link smearing with respect to an over-improved Symanzik 
type action ^| on the other. These two methods have their respective advantages compared to the 
approximate solution [2l|, [22[ of the Ginsparg- Wilson relation and APE smearing [28] . 

In this paper, the comparison has been made more complete and detailed, based first on the density- 
density two-point function and second on a point-by-point matching of the respective topological den- 
sities. The correspondence between the ultraviolet cutoff Xcut of the overlap analysis and the number 
of smearing steps justifies the use of over-improved stout-link smearing, which is computationally less 
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FIG. 16: The scalar (middle) and pseudoscalar (bottom) density of the first non-zero mode for the Q = — 1 
configuration shown previously in Fig. [131 along with the gluonic topological charge density after 200 sweeps of 
smearing (top). One sees how the scalar density extends over objects of differing charge, but that the regions 
of alternate charge are realized by the local chirality of the pseudoscalar density. In color online: negative 
density is blue/green, positive density is red/yellow (top), positive chirality is blue/green and negative chirality 
is red/yellow (bottom). In grey-scale: negative density dark, positive density light (top), and positive chirality 
dark, negative chirality light (bottom). 
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FIG. 17: The scalar (middle) and pseudoscalar (bottom) density of the first non-zero mode for the (5 = 1 
configuration shown previously in Fig. I14[ along with the gluonic topological charge density after 200 sweeps 
of smearing (top). Again, the scalar density extends over objects of differing charge, and the regions of 
alternate charge are revealed by the local chirality of the pseudoscalar density. In color online: negative 
density is blue/green, positive density is red/yellow (top), positive chirality is blue/green and negative chirality 
is red/yellow (bottom). In grey-scale: negative density dark, positive density light (top), and positive chirality 
dark, negative chirality light (bottom). 
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FIG. 18: The scalar (middle) and pseudoscalar (bottom) density of the first non-zero mode for the Q = — 1 
configuration shown previously in Fig. I15[ along with the gluonic topological charge density after 200 sweeps 
of smearing (top). Again, the scalar density extends over objects of differing charge, and the regions of 
alternate charge are revealed by the local chirality of the pseudoscalar density. In color online: negative 
density is blue/green, positive density is red/yellow (top), positive chirality is blue/green and negative chirality 
is red/yellow (bottom). In grey-scale: negative density dark, positive density light (top), and positive chirality 
dark, negative chirality light (bottom). 
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demanding, in the analysis of topological vacuum structure. This is of particular interest when inves- 
tigating the differences between the vacuum structure of quenched and non-quenched QCD [l5|, [2^. 

What we did not yet do in the present paper is the natural next step of formulating a description 
of the outstanding clusters of topological charge (at each level of ultraviolet smoothing) in terms 
of an filtered field strength, a possibility also offered by the overlap analysis [6]. Furthermore, one 
should study the degree of (anti)selfduality of the gauge field and the percolation properties of the 
corresponding (anti)selfdual domains [6]. Both questions are accessible not only by the fermionic 
method but also by the gluonic method based on stout-link smearing. 
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